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ABSTRACT 

AM Canum Venaticorum (AM CVn) binaries consist of a degenerate helium donor and a helium, 
C/O, or O/Ne WD accretor, with accretion rates of M — 10~ 13 — 10~ 5 Mq yr -1 . For accretion rates < 
10~ 6 Mq yr" 1 , the accreted helium ignites unstably, resulting in a helium flash. As the donor mass and 
M decrease, the ignition mass increases and eventually becom es larger than the donor mass, yielding 
a "last-flash" ignition mass of < 0.1 Mq. Bil dsten et al.l (|2007l ) predicted that the largest outbursts of 
these systems will lead to dynamical burning and thermonuclear supernovae. In this paper, we study 
the evolution of the He-burning shells in more detail. We calculate maximum achievable temperatures 
as well as the minimum envelope masses that achieve dynamical burning conditions, finding that AM 
CVn systems with accretors > 0.8 Mq will undergo dynamical burning. Triple-a reactions during 
the hydrostatic evolution set a lower limit to the 12 C mass fraction of 0.001 — 0.05 when dynamical 
burning occurs, but core dredge-up may yield 12 C, 16 0, and/or 20 Ne mass fractions of ~ 0.1. Accreted 
14 N will likely remain 14 N during the accretion and convective phases, but regardless of 14 N's fate, the 
neutron-to-proton ratio at the beginning of convection is fixed until the onset of dynamical burning. 
During explosive burning, the 14 N will undergo 14 N(a, 7) 18 F(a,p) 21 Ne, liberating a proton for the 
subsequent 12 C(p, 7) 13 N(a,p) 16 reaction, which bypasses the relatively slow a-capture onto 12 C. 
Future hydrodynamic simulations must include these isotopes, as the additional reactions will reduce 
the Zel'dovich-von Neumann-Doring (ZND) length, making the propagation of the detonation wave 
more likely. 

Subject headings: binaries: close — novae, cataclysmic variables — supernovae: general — white dwarfs 



1. INTRODUCTION 

Helium accretes onto helium, C/O, and O/Ne white 
dwarfs (WDs) at accretion rates ranging from M = 
10~ 13 — 10~ 5 Mq yr" 1 in AM Canum Venaticorum 
(AM CVn) binaries. iBildsten et all (|2007l ) showed that 
the resulting unstable helium flashes for M's between 
10~ 8 — 10~ 6 Mq yr -1 may lead to a new kind of ther- 
monuclear supernova (SN) that rises and falls in ~ 10 
days and reaches at most one-tenth the brightness of a 
normal Type la SN. These .la SNe have yet to be de- 
tected (the rec ent SN 2008ha com es the closest to resem- 
bling a .la SN: lFolev et al.1 12009') . but upcoming surveys 
such as the Palomar Transient Factory, Pan-STARRS-1, 
and Sky-Mapper will have cadenc es, depths, and sky co v- 
erage adequate to detect them if IBildsten et al.1 (|2007f )'s 
predicted rates in old stellar populations of ~ 1 per 10 4 
yr in a 10 11 Mq galaxy are correct. 

While it is known that the accumulated helium layer 
on the WD bur ns in an unstable ma nner once M < 
10~ 6 Mq yr" 1 (jlben fc Tutukovl llQSflh . it remains un- 
certain as to how the resulting flash evolves. In this 
paper, we calculate the time evolution of fully convec- 
tive He-burning envelopes as a series of hydrostatic so- 
lutions. This convective burning phase lasts for months 
— decades, much longer than the time it takes for trans- 
verse equilibration. Hence, we assume spherical sym- 
metry throughout our work. Due to hydrostatic radial 
expansion, the pressure at the base of the accumulated 



shell decreases, which leads to a maximum temperature 
for the base of the He-burning shell. For each WD mass, 
there is a minimum helium shell mass for which the nu- 
clear heating timescale becomes shorter than the time 
for a sound wave to traverse the radial thickness of the 
shell. Though this is not a definitive condition for the on- 
set of a rapid detonation, it certainly is the beginning of 
dynamical burning, either via a deflagration or a detona- 
tion. These minimum shell masses are na turally reached 
in th e evolution of the AM CVn binaries (Bildsten et al. 
2007). Envelopes smaller than these do not result in 
hydrodynamic burning as we describe above but instead 
complete their burning in a manner that allows for an ac- 
curate calculation with a hydrostatic code. Though mass 
loss occurs in these helium novae, like that which likely 
caused the outburst of V445 Puppis (iKato fc Hachisul 
120031 lliiima k Nakanishil 120081 IKato et al.ll2008[ ). such 
an event is not dynamically explosive. 

Our work provides the initial data relevant to the hy- 
drodynamic calculations needed to determine the out- 
come of these unstable flashes. In $2j we describe previ- 
ous calculations of helium ignition masses that motivate 
our range of envelope masses. In we describe our hy- 
drostatic calculation and calculate the maximum temper- 
atures reached during the envelope's hydrostatic evolu- 
tion, as well as the minimum envelope masses required to 
achieve dynamical burning. For AM CVn binaries with 
accretor masses > 0.8 Mq, the last (and largest) helium 
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flash will undergo dynamical burning. While the focus 
of this paper is on the evolution of He-burning shells in 
AM CVn binaries, our method's simplicity also lends it- 
self to other applications, such as the helium core flash 
at the tip of the red giant branch (RGB), which we ex- 
plore in §3.21 The effect of composition on our results is 
described in 21 We find that conditions for a-chain and 
14 N burning are not reached prior to dynamical burning, 
leaving unburned 12 C, 14 N, 16 O, and 20 Ne that will make 
a propagating detonation more likely in the thin shell as 
compared to calculations assuming pure helium. The nu- 
clear burning will be even further enhanced by the cat- 
alyzing effect of protons released in 18 F a-captures. We 
also find that the neutron-to-proton ratio is fixed between 
the beginning of convection and the onset of dynamical 
burning. In £J5J we summarize our work. 

2. HELIUM IGNITION MASSES ON ACCRETING 
WHITE DWARFS 

Though the purpose of this paper is to study the evo- 
lution of unstable helium flashes on helium, C/O, and 
O/Ne WDs, we must first set the stage by discussing 
the expected ignition masses, M\ m . All work done to 

date has been for constant M's and a range of core tem- 
peratures, T c , which are not necessarily the appropriate 
conditions for AM CVn binaries. However, these initial 
works prove adequate to reveal the parameters relevant 
for the AM CVn's. 

2.1. Numeric work and analytic expectation 

There have been many numerical calcula- 
tions of the ignition masses of H e- accreti ng 
WDs (INomoto fe Sugimotol \wfft iTaarf] fl980air 



Nomotol Il982at iFuiimoto fc Sugimotol |198S 



Woos lev. Taam. fe Weaverl 11 986; Ibcn & Tutukov 1991; 



Limongi & Tornambc 1991; Wooslcy & Weaver 1994: 



Piersanti et all 11 999; 



Piersanti. Cassisi. fc Tornambel 



20011; iGorielv et al.ll2002HYoon fc Langerll2004D . covering 
a large range of WD core masses (M c = 0.4 — 1.28 Mq) 
and accretion rates (M = 5x 10~ 10 -2x 10~ 7 M Q yr" 1 ). 
Note that the conductive opacity used in all of the pre- 
1991 stu dies was erro neou sly high due to a typogra phical 
error in llbenl (|1975f ). As llben fc Tutukovl (|1991f ) note, 
a higher conductive opacity reduces the heat flow into 
the degenerate core, so the envelope stays hotter than it 
should be, and thus the pre-1991 ignition masses are as 
much as a factor of 2 lower than they should be. 

Some authors only publish the accreted mass at ig- 
nition, Mj gn , which can be larger than the mass ini- 
tially contained in the convective zone due to the fact 
that the temperature peak in the hot accreting enve- 
lope can be above the envelope base. Thus, for M > 
3 x 10~ 8 Mq yr -1 , the co nvective envelope can be ini- 
tially as small as 0.5 Mj gn ([Fuiimoto fc Sugimotol [l~982t 
iLimongi fc Tornambel I1991D , and so M; gn is an upper 
limit to the convective envelope mass, M cnv . For M < 
10~ 8 Mq yr -1 , the envelope is not drastically hotter than 
the core, so the temperature peak at ignition is near the 
envelope base, and Mg nv ~ Mj gn (jLimongi fc Tornambel 
fT99lt llben fc Tutukovl [l99lh . 

Figure [T] displays numerically calculated ignition 
masses as a function of M, presumed to be constant, for 
0.6 Mq (squares), 0.8 M Q (triangles), and 1.08-1.2 M 




10" 8 
M [M Q yr" 

Fig. 1. — The points show M- lgn as a function of M for 0.6 Mra 
(squares), 0.8 Mq (triangles), and 1.08 — 1.2 Mq WDs (crosses) 
from numerical studies in the published literature; see text for 
sources. The values for the 1.08 — 1.2 Mq models may be too 
low due to an opacity error (see text). The donor mass evolution 
of a degenerate helium donor from C. Deloye (2007, private com- 
munication) is shown as a solid line; the donor mass evolves from 
upper-right to lower-left. Ignition masses expected for AM CVn bi- 
naries are within the shaded region. The evolution of a He-burning 
donor is shown as a dotted line, evolving first to higher and then 
to lower M jYungelsonl 120081') . 



(crosse s ) WD s. The 1.08 — 1.2 Mq, ignition masses 
dTaaml Il980bl; iNomotol Il982al: iFuimioto fc Sugimotol 
Il982t IWooslev. Taam. fc Weaverl Il986f ) were calculated 
before 1991, so they may be lower than an updated 
calculation would find. Note that the numerical igni- 
tion masses are nearly independent of M c for 0.6 Mq < 
M r < 1.2 Mq and M > 10~ 8 Mq yr" 1 . The 



data points shown in Figure [T] do not include the 
NCO chain ([Hashimoto et a l. 1986; P iersanti et~aL1ll999l; 
iPiersanti. Cassisi. fc Tornam bc 2001), which can ignite a 
flash prior to the onset of triple-a burning but is likely 
n egligible for AM CVn s ystems; see E|4. II for details. 

llben fc Tutukovl (fl989h used a one-zone calculation to 
estimate the helium ignition mass as 



M ign ~ 0.04 Mq 
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where Rb is the radius at the base of the accreted enve- 
lope. This rough calculation underestimates the ignition 
masses from the numerical work listed above by a factor 
of ~ 5 : or even a factor of 10 for M c > Mq. However, the 
scaling of M 1&n cx M~ 57 roughly matches the numerical 

results for M > 10~ 8 M yr -1 . 

In a ddition to their numerical work, llben fc Tutukovl 
(|1991f ) attempted to improve their previous one-zone an- 
alytic estimate by accounting for the inward luminosity 
into the core from the envelope. Their result is 
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where the helium ignition temperature in units of 10 s K, 
Tg,, is roughly equal to 1, and (3 is the fraction of the 
gravitational energy released by the accreting envelope 
material falling that goes into heating the core. Note that 
this M; gn is independent of M, which is in conflict with 

the numeric work listed above for M > 3x 10~ 9 Mq yr" 1 . 
However, the more gradual scaling with M c and Rb is in 
better agreem ent with the literature than the analytic 
estimate from llben fc Tutukovi (|1989fh 

2.2. AM CVn ignition masses 

The evolution of a fully-degenerate helium donor, 
as appropriate for AM CVn systems, is shown as a 
solid line in Figure [1] (C. Deloye 2007, private com- 
munication), beginning at M > 10~ 6 Mq yr" 1 and 
evolving to lower accretion rates. For M > 3 x 
10~ 6 M Q yr" 1 , outside of the bounds of Figure [TJ 
the He-burning in the accreted shell is thermally sta- 
ble (llben fc Tutukovi fp8l iTutukov fe Yungelsonl [l99l 
iShen fc Bildstenl 120071 ). potentially resulting in s hort- 
period supersoft sources (jvan Teeseling et~alTll997f) . As 
the donor loses mass, M decreases, the donor evolves 
from the upper-right of Figure [1] to the lower- left, an d 
the He-burning beco mes unstable (jBildsten et al.ll2006f ). 
iBildsten et alj (|2007l ) showed that the resulting ignition 
masses lie within the shaded region. When the donor 
mass becomes smaller than the required helium ignition 
mass, there will be no more helium explosions in the 
system. Thus, there is a "last-flash" M; gn that is the 
maxim um ignition mass for a given M c in an AM CVn 
system (IBildsten et al.ll2~007D . which Figure [1] shows to be 
< 0.1 Mq. Accurate M; gn values await calculations with 
M-evolution during the period of accumulation. These 
must also include the core heatin g from earlier phases of 
stable burning and weak flashes (|Tutukov fe Yungelsonl 
[1991 IBildsten et al.ll200l) . 

The AM CVn last-flash M ign is lower than the 
0.2 — 0.6 Mq helium shells expected in binary systems 
with He-burning donors, which were studied as potential 
Type la SN prog e nitors via the double-detonation sce- 
nario (iNomotoT Il982alfbl: IWooslev. Taam. fc Weaverl 
198(1 iLivnel Il99tt ILivne fc Glasnerl 119901 lT99lF 
Wooslev fc Weaver] 119941) . The donor evolution for 



a He-burning star i s shown as a dotted line in Figure 
Q] ()Yungel son 2008). The initial donor and accretor 
masses are 0.4 and 0.6 Mq. The system comes into 
contact at an orbital period of 80 min and evolves first 
to higher M's and then to lower M's under the influence 
of gravitational wave radiation. The M-evolution for 
these stars leads to Mj gn > 0.2 Mq, significantly higher 
than those expected in AM CVn systems. The strong 
overlap of the calculated points in Figure [1] with the 
He-burning star evolution is due to the substantial work 
on double-detonations. 

3. CALCULATION AND EXAMPLE 

AM CVn ignition masses can be as large as 0.1 Mq, 
and because the accreted envelope can be such a large 
proportion of the WD mass, it is important that 



the envelope and core variables are calculated consis- 
tently so as to maintain hydrostatic equilibrium. Our 
calculations thus differ fro m the analytical work of 
iFuiimoto fc Sugimotol (|1982f ). who assumed that the ra- 
dius of the underlying WD remains unchanged through- 
out the evolution of the flash. 

The lowest density in the underlying WD, just un- 
der the base of the accreted envelope, is ~ 10 5 g cm" 3 
at the time of helium ignition. The associated Fermi 
energy is equal to the thermal energy when T — 4 x 
10 8 K (p/10 5 g cm" 3 ) 2 / 3 . Thus, for a core temper- 
ature T c < 10 8 K, the temperature will not signifi- 
cantly affect the core structure. For simplicity, we choose 
T c = 10 7 K for all of our models (unless otherwise 
noted as in i)3.2|) as expected for the majority of the 
helium flashes and a factor of a few lower than the max- 
imum 2^_e2crje£ted_for_the last flash in AM CVn sys- 
tems (jBildsten et al.ll2006l) . The core is assumed to have 
a homogenous composition of 50% C and 50% O by 
mass, but the core structure is nearly independent of 
the composition choice because 4 He, 12 C, 16 0, and 20 Ne 
all have the same charge-to-mass ratio. The equation 
of state (ISaumon et aTTl 19951 : iTimmes fc Swestvl 120001: 
Rogers fc Navfonovl l2002f) " opacity (llglesias fc Rogers! 
19931 I1996D . nuclear burning network (iTimmesI I1999L 
unles s otherwise noted), neutri no cooling (lltoh et alj 
19961). and electron screening (|Gr aboskc e t al.l 11973 



Alastuev fc Jancovicil Il978t lltoh et al.lll979f ) are calcu- 



lated with the MESA code package. 1 

The density and radius are integrated outwards in mass 
coordinates from a central density, p c , under the assump- 
tion of hydrostatic equilibrium until the desired WD core 
mass, M c , is reached. Specifying p c , M c , and T c gives the 
pressure and radius at the WD core-envelope interface. 
Since the pressure and radius are continuous variables, 
these are equal to the pressure and radius at the base of 
the accreted envelope, P, and Rb, whereas the tempera- 
ture and composition are discontinuous at the WD core- 
envelope boundary by construction. The assumption of 
a discontinuous composition relies on minimal diffusion 
during the accretion phase, justified in Appendix 

The actively-burning envelope is an efficient convec- 
tion zone, 2 so that the temperature profile is given by 
the adiabatic index, V a d = <91nT/<91nP| s . The pressure 
and radius at the base of the envelope, calculated from 
the core integration, are used as the starting point of the 
envelope integration. A base temperature, is speci- 
fied, and the envelope variables are integrated outwards 
until P(r) <C Pb, yielding the envelope mass, M onv , as a 
function of Tb for a particular combination of Pb and Rb. 
For a given M c , there is at most one Tf, for each Pf, that 
yields the desired M cnv . Thus, for fixed M c and M env , 
we can describe the evolution of Pf, and Tb during the 
flash as the envelope's entropy is increased by nuclear 
burning. 

Figure [2] shows the evolution of the base temperature 
and pressure as the entropy increases for 0.01, 0.1, and 
0.2 Mq envelopes on a WD with M c = 0.6 Mq (solid 
lines). Figure [3] shows the same, but for 0.01, 0.05, 

1 http://mesa.sourceforge.net/ 

2 The convective speed, v a> can be defined via the convective flux 
as F c ~ pv'c . Efficient convection is equivalent to the condition that 
v c is much smaller than the sound speed. 
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Fig. 2. — Hydrostatic evolution of T b and Pj, for fully convective 
0.01, 0.1, and 0.2 Mq pure helium envelopes on a 0.6 Mq WD 
with T c = 10 7 K (solid lines). The base temperatures at ignition 
are ~ 10 s K. Conditions for which the global heating timescale, 
t heat , equals various timescales up to 1 minute are labelled. As 
the envelopes evolve to higher entropy configurations, T b rises at 
nearly constant pressure until the scale height becomes ~ 20% of 
the WD radius, at which point the envelope expands and Pf, and 
T;, decrease. The radiation pressure limit, P b = (dotted 
line), provides an upper bound to T b . Also shown are conditions 
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Fig. 3.— Same as Fig. [2] but for 0.01, 0.05, and 0.1 Mq envelopes 
on a 1.0 M e WD. 



and 0.1 Mq envelopes on a 1.0 M Q WD. The envelopes 
are composed of pure helium for simplicity because the 
triple-a reaction is the dominant entropy source for dis- 
playing the envelope evolution and dynamical burning 
conditions that are shown in Figures [2] and [3l see <j4] for 
details on the effects of envelope composition. 

The hydrostatic growth of the envelope is governed by 
a global heating timescale ihcat = J (cpTdm) / J(edm), 
where cp is the specific heat, and the energy generation 
rate is e. This timescale is not to be confused with the 
local heating timescale t' heat = cpTb/e; see jj3.il The in- 
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Fig. 4. — Maximum T b for convective helium envelopes. Con- 
tours of constant maximum T b are labelled by their logarithmic 
value. The minimum M cnv that achieves dynamical burning (as 
defined by reaching t^ yn = tC eat ; see £13-16 is shown by the dotted 
line. Envelopes larger than this minimum M env become dynamical 
before reaching their maximum T b . Lowering the requirement to 

shown in Fie 



0.1t[ lcai _reduces the minimum M env by a factor of ~ 2, 



tegrals are performed over the entire convective envelope 
because even though the majority of nuclear burning oc- 
curs at the base of the envelope, th e entropy release must 
heat the whole con vective zone (jWeinberg et al.l 120061 : 
iPiro fc Chanel [2~008f ) . Conditions for which ih eat equals 
a year, a month, a day, an hour, and a minute are la- 
belled in Figures and [31 if the envelope evolves quickly 
enough. Note that iheat first decreases with increased en- 
tropy but then increases after radial expansion becomes 
important and Tf, reaches its peak. 

Convection is initiated when Tt ~ 10 8 K. As He- 
burning increases the entropy of the envelope, T b rises at 
nearly constant pressure because the envelope is initially 
geometrically thin; the ratio of pressure scale height, h, 
to core radius, assuming pressure due to non-relativistic 
electron degeneracy, is 



h ™ 

Kb 



Pb 



10 22 dyne cm 2 
M Q \ ( R b 



2/5 



5 x 10 8 



(3) 



where h = Pb/pbg and g = GM c jR\ is the gravitational 
acceleration at the envelope base. Thus, Pb remains 
nearly constant during the initial phase of the flash and 
is approximately equal to GM c M env / AnRf. 

The ideal gas pressure and radiation pressure become 
important as the entropy and temperature increase, and 
the envelope becomes geometrically thick, resulting in a 
decrease of Pb- Further entropy increase when the shell 
is thick results in a decrease of Tb, with the maximum 
temperature achieved when h ~ 0.2Rb- This maximum 
Tb only depends on M c and M cnv for a given envelope 
composition. Figure [3] shows contours of maximum Tb for 
pure helium envelopes as a function of the WD mass. The 
contours are truncated at the M env 's above which the 
burning becomes dynamical before the maximum tern- 
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perature is reached (dotted line); see ij3.ll 

Because of the ultracompact nature of AM CVn bina- 
ries, we must consider the possibility of binary interac- 
tions during the evolution of the helium flash. At or- 
bital periods > 3 min, as expected during the helium 
flash regime when M = 3 x 10~ 8 - 2 x 10 ~ 6 Mq yr _1 
(|Nelemans et all 120011: iBildsten et aLll2006f) . the binary 
has a separation > 4 x 10 9 cm, wide enough that the 
radius of the burning WD stays within its Roche radius 
for envelopes that achieve dynamical burning. However, 
the evolution of non-dynamical envelopes after the peak 
temperature must account for the impact of Roche-lobe 
filling and a possible common envelope phase. We leave 
these details for future work. 

3.1. Dynamical burning 

Our calculations have assumed that the convective mo- 
tion is very subsonic (i.e., the convection is efficient), so 
that the envelope is in hydrostatic equilibrium. How- 
ever, if the base of the envelope heats up more rapidly 
than it can adjust to maintain hydrostatic equilibrium, 
this assumption is contradicted. To quantify this condi- 
tion, we compare the local heating timescale at the base, 
^hcat = c pTb/tb, to the dynamical timescale, idyn = h/c s , 
where the pressure scale height, the sound speed, c s = 

dP/dp\ ad , the specific heat, and the energy genera- 
tion rate are evaluated at the envelope base. The local 
timescale, t' heat , as opposed to the global timescale, ihcat, 
is the appropriate choice because we wish to know when 
the base conditions decouple from the rest of the enve- 
lope. The majority of nuclear burning occurs near the 
base of the envelope, but the thermal energy content is 
spread out in the envelope, so t' hcat is shorter than ih ca t 
by a factor of ~ 10. 

When the local heating timescale becomes shorter than 
the dynamical timescale and the envelope cannot ad- 
just itself to maintain hydrostatic equilibrium, an over- 
pressure devel o ps, w hich may lead to a detonation. As 
IBildsten et alj (|2007t ) explain, such a detonation at the 
relatively low ignition pressures achieved in AM CVn sys- 
tems produces radioactive isotopes such as 48 Cr, 52 Fe, 
and 56 Ni. The combination of their short half-lives and 
the small ejecta masses yields explosions with bright- 
nesses comparable to sub-luminous Type la SNe, but 
with more rapid decays of ~ 5 — 10 days. 

Figures [2] and [3] show conditions for which tdyn = t' heat 
for a 0.6 and a 1.0 M© WD, respectively (dashed lines). 
We also show conditions for which td yn — O-l^deat m or der 
to account for uncertainty in the exact requirements for 
the initiation of dynamical burning. Since dynamical 
burning is achieved before an envelope reaches its peak 
Tb, the contours in Figure 0] are truncated for envelopes 
that reach t dyn = t' hoat . 

It is clear that large envelopes will burn dynamically, 
while smaller ones will not, and that there is a mini- 
mum envelope mass for which dynamical burning con- 
ditions are reached. Typical timescales at the onset of 
dynamical burning for these minimum mass envelopes 
range from idyn = i nca t = 0.1 — 0.5 s. Figure [5] shows 
these minimum dynamically-burning M env 's vs. WD 
mass for bot h td y n = ^heat ano - (solid lines). As 

discussed in £12.21 the values of the "last-flash" Afi gn 's 
for AM CVn systems remain somewhat uncertain due 
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Fig. 5. — Minimum M en v's that reach dynamical burning con- 
ditions vs. M c (solid lines). To account for uncertainties, masses 

that achieve tdyn = * nca t an< ^ *dyn = 0-l*hcat are snown - The 
shaded region shows M cnv = 0.03 — 0.1 Mq, which is an es- 
timate of the "last- flash" M- lgn expected for AM CVn binaries 
UBildsten et al.ll2007h . 

to the as-yet uncalculated effects of binary evolution 
and core heating. Thus, we show a range of last-flash 
M[ sn = 0.03 — 0.1 Mq (shaded region), which indicates 
that accretor masses > 0.8 Mq will undergo dynami- 
cal burning in AM CVn systems, giving rise to .la SNe. 
Large WD masses also likely demand O/Ne WDs, which 
would obviate any concern of a carbon detonatio n deep 
in the accreting WD t riggered via shock focusing (|Livnd 
[19901 : iFink et~aLl[2007h . 

3.2. Helium core flash on the red giant branch 

The simplicity of our calculation can be applied to any 
situation involving a convective He-burning shell on an 
isothermal core, with the caveat that the outer boundary 
condition at the edge of the convective shell may be dif- 
ferent for different circumstances. In particular, we now 
briefly examine the helium core flash at the tip of the red 
giant branch (RGB), which produces convective helium 
burning shells on a partially degenerate helium core. 

As a < 2Mq star evolves up the RGB, a tem- 
perature inve rsion develops i n the growing degenerat e 
helium core (jThomad Il967t ISweigart fc Gross! Il978h . 
When the helium core m ass reaches 0.48 Mq (see 
iSerenelli fc Fukugital 120051 for recent calculations) , the 
temperature peaks a t a mass coordinate of 0.3 Mq 
(jFuiimoto et al.lfT990l ). yielding a convective He-burning 
shell with M c = 0.3 Mq and M env = 0.18 Mq (e.g., 
ISweigart fc Gross! 119781 : llben fc Renzlnl 1 1984ft. T his is 
the first of ~ 10 core flashes jlben fc Renzinil ll984). with 
the mass coordinate of the convective envelope base mov- 
ing inwards with each flash. In ~ 10 6 years, the degener- 
acy has been removed from the helium core, and the star 
undergo es stable helium core b urning on the horizontal 
branch (|Serenelli fc Weisdl2005ft . 

Figure [6] shows contours of peak temperature as in Fig- 
ure |4j but for helium envelopes on a hot helium core with 
T c = 10 8 K as appropriate for these helium core flashes. 
While the exact core composition and T c do not matter 
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Fig. 6. — Same as Fig. [4] but for helium envelopes on helium 
cores with T c = 10 s K as appropriate for helium core flashes at the 
tip of the red giant branch, and with peak temperature contours 
ranging from 10 8 ' 2 — 10 8 ' 7 . The dashed line shows a constant total 
mass of M c + M env = 0.48 Mq , and the cross shows the first core 
flash, which results in a 0.18 Mq convective envelope on a 0.3 Mq 
helium core. 

for our other calculations, they are important for these 
low mass calculations with M c < 0.4 Mq because the 
small hot core is only partially degenerate. The dashed 
line shows values of M c and M onv that give a constant 
helium mass of M c + M cnv = 0.48 Mq, the mass of the 
helium region at the tip of the RGB, and the cross shows 
the core and envelope mass of the first helium core flash. 
Our calculation does not take into account the relatively 
high entropy boundary condition at the top of the he- 
lium convective zone, which results from the hydrogen- 
burning shell above. However, this boundary condition 
has very little effect on the condit ions at the base of the 
helium layer; see, e.g., Figure 8 of lThomasI (fl967h . 

Figure [5] shows clearly that the peak temperature of all 
of the flashes is ~ 2 x 10 8 K. Moreover, we see that these 
flashes never become dynamical, since the combination of 
M c and M cnv yielding a total mass of 0.4 8 Af B is always 
in th e hydrostatically-burning regime (jFujimoto et al.l 
11990ft . Thus, helium" core flashes should be accurately 
calculable by presuming hydrostatic balance and effi- 
cient convection. In fact, summing M c and the mini- 
mum dynamically-burning M env in Figure [6] shows that 
no pure helium star < 0.6 Mq can ever achieve dynami- 
cal burning. This calculation clarifies wh y modern multi- 
dimensional hydrodynami c calculations (jDearborn et alJ 
120061: iMocak et alll2008al |b1) nearly always find that the 
flash leads to convective velocities close to that expected 
from mixing length and temperature gradients close to 
the adiabatic gradient. 

4. ENVELOPE COMPOSITION 

The previous calculation assumed that the convective 
envelope consists of pure helium, neglecting the pres- 
ence of other elements that have been accreted from the 
donor, produced by nucleosynthesis, or dredged up from 
the underlying WD core. These trace elements do not 
greatly affect the Pb — Tf,-trajectories of the envelopes be- 
cause those only depend on the equation of state, which 



is always dominated by helium. The composition also 
does not strongly affect the heating timescales, as the 
trace elements are depleted when burning conditions are 
reached. 

However, the propagation of detonations may be af- 
fected by the presence of 14 N, 16 0, 18 F, 20 Ne, and 
24 Mg, which lead to nuclear reactions that are more 
volatile than the triple-a process at post-shock temper- 
atures > 10 9 K. These isotopes can accelerate the burn- 
ing and shorten the Zel'dovich-von Neumann-Doring 
(ZND) detonation length. All pr evious work has only 
been for pure helium comp osition (jKhokhlovil 19881 119891 : 
iTimmes fc Niemeverl l2000i ) and thus neglects these ef- 
fects. Furthermore, the occurrence of weak interactions 
changes the neutron-to-proton fraction, impacting the 
nucleosynthetic outcome of the detonation. In this sec- 
tion, we explore the range of abundances expected at the 
onset of dynamical burning in order to motivate future 
hydrodynamical calculations. 

4.1. Accreted composition 

The material that is accreted from the degenerate he- 
lium donor in an AM CVn system is nearly pure helium. 
Assuming the donor has undergone CNO-burning of its 
hydrogen, the majority of the accreted CNO isotopes will 
exist as 14 N because the 14 N proton-capture is the slow- 
est step of the CNO cycle for typical CNO-burning con- 
ditions. Thus, we expect the accreted envelope to have a 
14 N mass fraction of Xu n — 9 x 10~ 3 ( the sum of solar 
CNO mass fractions from lLodde rs 2003), a helium mass 
fraction of Y = 0.99, and a negligible amount of heavier 
elements. 

Because the triple-a process is active at lower tem- 
peratures than a-captures on 14 N, the mass frac- 
tion of 14 N will not change prior to the onset 
of thermally-unstable triple-a burning and convec- 
tion unless conditions are reached for 14 N electron- 
captures, which are the beginning of the NCO 
chain: 14 N(e~, i/) 14 C(a, 7 ) 18 0(a, 7) 22 Ne. However, 
the threshold density required for the e~-captures to 
outpace the reverse /3-decays is 1.2 x 10 6 g cm~ 3 
for a He-dominated plasma ([Hashimoto et all 119861 : 
iPiersanti. Cassisi. fc Tornambel 120011 ). which corre- 
sponds to an electron degeneracy pressure of 4.0 x 
10 22 dyne cm~ 2 . For 0.65 Mq < M c < 1.3 M , the 
minimum M env 's that achieve idyn = t'h ca t begin their 
convective evolution with Pb lower than this threshold 
pressure, and thus A14N will remain at 10~ 2 . For en- 
velopes with Pb > 4.0 x 10 22 dyne cm -2 at the onset of 
convection, it is possible that some of the accreted 14 N 
will capture electrons. However, the timescales associ- 
ated with the 14 N electron-capture are > 10 7 yr, and 
thus it is unclear if there is enough time for the conver- 
sion of 14 N to 14 C in the 10 6 yr it takes to accrete 0.1 M Q 
at 10 -7 Mq yr -1 even if the threshold density is reached. 
Furthermore, if electron-captures do occur, they will oc- 
cur only at the bottom of the accreted envelope, where 
the density is highest. A proper calculation of the 14 N 
and 14 C abundances awaits a time-dependent simulation 
and thus we explore the fate of 14 C in i)4.3.2[ but it is 
likely that the majority of accreted 14 N nuclei will still 
be 14 N at the onset of convection in an AM CVn system. 

4.2. Composition changes due to core dredge-up 
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Analytic, nu meric, and observational work on classi- 
cal novae (e.g.. IStarrfield et al.lll972t iPrialnik fe Kovetd 
fl98llTruran fc Liviolll986HGehrz et al.lll998l ) show that 
core material is often entrained in the convective enve- 
lope. The amount of core material can be as high as 
30% by mass in classical nova ejecta, and such a dredge- 
up would have a large effect on the abundances of 12 C, 
16 0, or 20 Ne in the convective helium envelope, depend- 
ing on the composition of the C/O or O/Ne core material 
directly underneath the accreted envelope. 

However, the exact process by which the classical nova 
env elope becomes enriched in core material is uncertain 
fsee lShen fc Bil dstcn 2009 for a summary), and it is not 
clear that the same mechanisms also function for helium 
flashes. For example, one possible classical nova core- 
envelope mixing process invol ves the diffusion of hydro - 
gen into 12 C-rich core material (|Prialnik fc Kov ctz 1984). 
Convection is then initiated below the bulk of the ac- 
creted layer in the diffusing hydrogen tail via p+ 12 C re- 
actions. However, in the case of helium accretion, the 
diffusion of helium into 12 C-, 16 0-, or 20 Ne-rich material 
below the layer will not trigger a flash because the re- 
action rates for a-captures on these elements at ~ I0 8 
K are much slower than the triple-a reaction rate. Con- 
vection is thus almost certainly initiated at or above the 
base of the accreted layer via triple-a or NCO reactions. 

Core-envelope mixing in classical novae is also hy- 
pothe sized to occur via convective overshoot (|Wooslevl 
119861) or shea r mix i ng du e to the convective motion 
(|Glasner et all H997L [20071 IRosner et all l200l . Such 
mixing mechanisms have not been studied in the con- 
text of convective helium shells, so it is unclear whether 
or not they can significantly alter the envelope composi- 
tion. Since core-envelope mixing is not decisively ruled 
out, we allow for the possibility of core-envelope mixing 
in gXUby considering the fate of 12 C, 16 0, and 20 Ne 
in the convective envelope. 

4.3. Compositional evolution due to nucleosynthesis 

We now calculate the evolution of 14 N, 14 C, 12 C, 16 0, 
and 20 Ne due to nuclear burning. While our methodol- 
ogy prohibits the time-dependent calculation of an iso- 
tope's abundance, we can determine its absence or pres- 
ence by comparing ihoat to its depletion timescale, ti = 
\d\nXifdt\~ 1 . The depletion timescale in a uniformly- 
mixed convective envelope with respect to a-captures is 



' J nn{av)i dm 

where the integral is performed over the convective en- 
velope, 77,4 is the number density of He, and n4{av}i 
is isotope i's a-capture rate. When an isotope is con- 
sumed more rapidly but is produced more slowly than 
the envelope evolves, its abundance drops exponentially. 
For example, if the lifetime of 16 O with respect to a- 
captures becomes shorter than ihcat but the lifetime of 
12 C remains longer than iheat, any 16 O present will be 
converted to 20 Ne (and then possibly to an isotope higher 
in the a-chain) without being replenished by the conver- 
sion of 12 C to 16 0. 

4.3.1. The fate of a-chain elements 



We begin by calculating the lower limit to the 12 C mass 
fraction, X\2, set by He-burning. Since the vast major- 
ity of heat is generated by the conversion of 4 He to 12 C, 
we can estimate the 12 C mass fraction, X±2, produced 
through nucleosynthesis at a given evolutionary stage by 
calculating the difference in total energy between that 
stage and that of the envelope at the onset of convection. 
This will be a lower limit to the actual 12 C abundance 
at that time because core-envelope mixing may intro- 
duce additional carbon from the core. Furthermore, we 
neglect the kinetic energy of convective motion in this es- 
timate, so the actual difference in energies and thus X12 
may be much larger. 

We must take care to include the energy contained in 
the core as well as the envelope because the significant 
radial expansion of the envelope during the course of the 
convective evolution means work is done on the envelope 
by the core. The total energy in the WD, neglecting the 
kinetic energy of convection, is 

-Etotai = J (u- —-rj dm -> ( 5 ) 

where u is the specific internal energy. The minimum 
Monv that achieves dynamical burning for a 0.6 Mq 
WD is 0.19 Mq. To evolve from an initial state of 
Tb = 10 s K to the stage where td yn = t'hc&t requires 
A-Etotal — 8 x 10 48 erg, which is made available by burn- 
ing 7 x 10~ 3 Mq of 4 He to 12 C, giving a final carbon 
mass fraction Xu — 0.04. For a 0.05 Mq envelope on a 
1.0 Mq WD, a similar analysis yields X12 ~ 0.05 at the 
onset of dynamical burning. For larger envelopes than 
these, the energy required to achieve dynamical burning 
is smaller (because Xfc at the onset of dynamical burning 
is lower; see Figs. [2] and [3]) and M cnv is larger, yield- 
ing X12 = 10~ 3 — 10~ 2 at the onset of dynamical burn- 
ing. However, as previously mentioned, these values are 
a lower limit due to our neglect of the material dredged 
up from the core and the energy in convective motions. 

We now consider the fate of trace a-chain elements 
heavier than helium (i.e., 12 C, 16 0, 20 Ne, and 24 Mg). 
The lines in Figure [7] show the base temperatures and 
pressures for which ti of equation 3] is equal to theat- The 
global isotope depletion and envelope heating timescales 
are functions of M c , but due to the strong temperature- 
dependence of nuclear burning, the isotopic depletion 
lines are nearly independent of M c ; for simplicity, we only 
show lines calculated for M c = 0.6 Mq. The unscreened 
a-capture lifetimes for 16 (dashed- dotted line) and 24 Mg 
(long-dashed line) are calculated with the MESA code, 
while the lifetime of 12 C (short-dashed line) is calculated 
from iKunz et all (|2002fl . The depletion line for 20 Ne is 
not shown but is nearly identical to that of 16 O. The 
radiation-pressure limit is represented by a dotted line. 
Also shown are the trajectories of 0.1 and 0.2 Mq pure 
helium envelopes on a 0.6 Mq WD (thick solid line). The 
trajectory of the 0.2 Mq envelope has been truncated 
at the dynamical burning line where idyn = ^hcat (thick 
dashed line). As with the depletion lines, the conditions 
for which idyn = t' hcat are a function of M c but depend 
only slightly on it because nuclear-burning is strongly 
temperature-sensitive. 

Isotope lifetimes are longer than iheat in the regions 
below the depletion lines, so an isotope's abundance does 
not decrease due to a-captures significantly in this area 
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Fig. 7. — Base temperatures and pressures for which isotope life- 
times, ti, are equal to theat- The isotopic lifetimes for 12 C (short- 
dashed line), le O (dashed- dotted line), and 24 Mg (long-dashed 
line) are calculated with respect to a-captures. The lifetimes are 
longer than theat f° r temperatures and pressures below the deple- 
tion lines. Also shown is the radiation pressure limit (dotted line) 
and base trajectories for 0.1 and 0.2 Mq pure helium envelopes on 
a 0.6 Mq WD (thick solid line), the latter of which is truncated at 
the dynamical burning line where id yn = *heat (*' l * c ' c dashed line). 

of parameter space, although it may increase if a different 
isotope is converted into it. In particular, the 12 C mass 
fraction does not decrease from a-captures because the 
hydrostatic evolution of all envelopes remains below the 
12 C depletion line, but X 12 does increase due to triple-a 
reactions. 

Figure [7] shows that depletion of a-chain isotopes 
occurs in a region of parameter space inaccessible to 
hydrostatically-burning convective envelopes, and thus 
they are still unburned when the envelope becomes dy- 
namical. If a shock does develop and the post-shock 
temperature is > 10 9 K, the unburned a-chain elements, 
and especially 16 0, 20 Ne, and 24 Mg, capture a-particles 
faster than the triple-a process operates. Thus, hydro- 
dynamic simulations of helium detonations must include 
unburned a-chain elements resulting from dredge-up and 
nucleosynthesis, as they can reduce the ZND length as 
compared to pure helium detonations. 

4.3.2. The fate of 14 N and 14 C in envelopes that become 

dynamical 

We now calculate the evolution of the accreted 14 N, 
a small amount of which will be converted to 14 C if 
Pb > 10 6 g cm -3 and the accretion timescale > 10 7 yr. 
Due to the /3-decays in their reaction chains, these iso- 
topes determine the neutron-to-proton ratio at the onset 
of dynamical burning, which will affect the isotopic yield 
of the detonation. Furthermore, as we will see, they can 
affect the progress of the helium detonation by providing 
protons to catalyze other reactions. These isotopes will 
also play a role in the nucleosynthetic outcome of en- 
velopes that do not achieve dynamical burning, but we 
leave these details for future work. 

Conditions for which ihcat equals the positron-emission 
lifetime of 18 F (dashed- dotted line) and the unscreened a- 
capture lifetimes for 14 N (long-dashed line), 14 C (short- 
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Fig. 8. — Same as Fig. [3 but with different isotopes: the isotopic 
lifetimes for 14 N (long-dashed line) and 14 C (short-dashed line) are 
calculated with respect to a-captures, while the lifetime of 18 F is 
shown with respect to both «-captures (solid line) and /3-decays 
(dashed- dotted line). 

dashed line, iGorres et al.l Il992h. and 18 F, which y ields 
21 Ne and a proton (solid line; iKarakas et aLll2008D . are 
shown in Figure [H] The other lines are as in Fig- 
ure [71 with the same caveat that while the calculations 
shown arc for M c = 0.6 Mq, they are nearly identical 
for other core masses due to the strong temperature- 
sensitivity of nuclear burning. It is clear that any 14 C 
that is produced prior to convection will remain as 14 C 
at the onset of dynamical burning. Decays back to 14 N 
have a half-life of 5730 yr (or even longer due to ex- 
clusion principle inhibition at the base of the en velope ; 
iPeterson & Bahcall|[l96l iHashimoto et al.lll986D , which 
is much longer than the convective burning phase. Fur- 
thermore, a-captures are negligible because the 14 C de- 
pletion line crosses the dynamical burning line at a pres- 
sure < 4 x 10 22 dyne cm -2 ; i.e., a-captures on 14 C at the 
high densities required to create 18 O occur at tempera- 
tures > 5.5 x 10 8 K, higher than necessary for dynamical 
burning at these densities. 

Likewise, 14 N will remain unchanged prior to dynami- 
cal burning in nearly all cases. If conditions are reached 
for 14 N(a, 7) 18 F reactions, the 18 F will not have a chance 
to /3-decay because its /9-decay half- life of 110 min is 
much longer than the timescale at the onset of dynam- 
ical burning, tdyn = ^hcat- Electron-captures onto 18 F 
during the hydrostatic envelope evolution can also be ne- 
glected, as their timescale is much longer than thsat. for 
cond i tions after 14 N a-captures hav e taken place (|Bahcalll 
[1961: [LaHgl[l980l iTillev et al.lll995h . Therefore, no weak 
interactions will occur prior to dynamical burning for nu- 
clei that are 14 N at the onset of convection, even if they 
are converted to 18 F. 

Thus, the neutron-to-proton ratio at the onset of dy- 
namical burning is determined at the beginning of con- 
vection: nucleosynthesis during the convective phase will 
not change this ratio before dynamical burning occurs. If 
14 C exists prior to convection, it will remain 14 C, and 14 N 
will either remain unburned or capture an a-nucleus to 
become 18 F, but will not undergo any weak interactions. 
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As with the a-chain isotopes, it is important to 
consider the presence of 14 C and 18 F when cal- 
culating the hydrodynamic evolution of the detona- 
tion. At temperatures > 10 9 K, the a-capture chain 
14 C(a,7) 18 0(a, 7) 22 Ne will take place, although its en- 
ergy input will likely be negligible due to the small 
abundance of 14 C. More significantly, a-captures on 
18 F lead to the production of protons via 18 F(a,p) 21 Ne 
(recently examined in the context of AGB evolu- 
tion in iKarakas et al.l [2008), which act cLS ct Celt - 
alyst for a-chain rea ctions via 12 C(p, 7) 13 N(a,p) 16 
([Weinberg et al.|[2006h . This process bypasses the com- 
paratively slow 12 C a-capture and will certainly impact 
the evolution of the detonation. Furthermore, this reac- 
tion will play a role in determining the nucleosynthetic 
outcome of envelopes that do not become dynamical, the 
details of which we leave to future work. 

5. CONCLUSION 

Motivated by iBildsten et alj (|2007l )'s result that dy- 
namical helium burning is naturally attained during the 
evolution of AM CVn binaries, we have studied the 
progress of convective helium envelopes via hydrostatic 
calculations under simplifying assumptions O- We 
have found the maximum possible temperatures (Fig. [5]) 
as well as the minimum envelope masses that achieve 
dynamical burning, as defined when the local heating 
timescale, t' heat , equals the sound-crossing time of the 
envelope, td yn (Fig. [5]). The generality of our calcula- 
tion allowed us to also examine the helium core flash at 
the tip of the red giant branch ( i)3.2p . confirming pre- 
vious results that the flash is a slow enough event that 



hydrostatic evolution is a safe assumption. 

In order to motivate the starting conditions for future 
hydrodynamic studies, we examined the nucleosynthetic 
fate of elements in <|4] We found that a-chain elements 
such as 12 C, 16 0, 20 Ne, and 24 Mg, produced via nucle- 
osynthesis or dredged up from the degenerate core, re- 
main unburned during the hydrostatic envelope evolution 
(Fig. [7]) and will be present if a detonation commences, 
affecting the ZND length. 

For large envelope and core masses, we have found 
that some 14 N can capture electrons during the accretion 
phase to form 14 C, but this is likely a small effect. In- 
dependently of envelope and core masses, we have found 
that the neutron-to-proton ratio at the onset of dynam- 
ical burning, which will determine the detonation's iso- 
topic distribution, is the same as at the onset of convec- 
tion because there is no opportunity for weak interac- 
tions during the convective evolution to dynamical burn- 
ing; the timescale for the obvious candidate (positron- 
emission by 18 F) is so long that envelopes that achieve 
dynamical burning do so before 18 F can /3-decay (Fig. 
[8]). The existence of 18 F is very important for determin- 
ing the ZND length because of the production of protons 
via 18 F(a, p) 21 Ne, which will bypass the slow step of the 
a-chain by catalyzing the conversion of 12 C to 16 O. 
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APPENDIX 

A. CHEMICAL DIFFUSION BETWEEN THE WD CORE AND ENVELOPE 



In this Appendix, we justify our assumption of a discontinuous composition profile at the WD core-envelope interface. 
Due to electron degeneracy, there is a negligible differentiating force on helium in a degenerate C/O or O/Ne plasma. 
Hence, the accreted helium simply diffuses into the WD core at the rate given by a random walk and penetrates 
a distance rj n ~ V Dt^ cc during the time it takes to accumulate an ignition mass, i acc = M- 1&n /M, where D is the 
diffusion coefficient. For a carbon plasma, the Coulomb coupling parameter of the one component plasma (OCP), which 

measures the strength of ion-ion interactions, is Pi2,i2 = l-7p^ 3 /T$, where p§ = p/10 5 g cm" 3 and Ts = T/10 8 K. 
Thus, the core-envelope interface is in an intermediate regime between weak and strong ion coupling, for which an exact 
value for D at these temperatures and densities has not yet been calculated for the case of helium diffusion in a carbon 

plasma ([Paquette et al.lll986[ K The fitting formula of llben fc MacDonaldl (|1985h yields D = 0.03 cm 2 s" 1 T 8 5/2 /p5, 
neglecting the temperature- and density-dependence of the Coulomb logarithm. This estimate yields 

T 5/4 , , -.1/2 

T~° m TjVfe) ■ < A1 > 

where the pressure scale height near the core-envelope interface is h = P&/ pbg, and g% is the gravitational acceleration, 
g, in units of 10 8 cm s~ 2 . Thus, diffusion is unimportant for t acc < 10 s yr, or M > 10 -9 Mq yr" 1 for Mj gn < 0.1 M r, 
as long as the diffusion coefficient is not actually orders of magnitude higher than given bv llben & MacDo nald (198j|'s 
fitting formula. 

An alternate estimate for D can be made by appealing to results for the OCP. No appropriate calculations are 
available for trace helium in a C/O OCP, but we could attempt to follow B ildsten a nd Hall (2001) by using a Stokes- 
Einstcin-like argument, assigning a finite size to the helium ion. However. IDaligaultl (2006) showed that the expected 
Stokes-Einstcin relation between Docp and the viscosity of the medium, rj, only begins to hold in the OCP once 
F > 40. Hence, our only real choice for an alternate estimate is to use an approximate result relevant for T ~ 1 — 10, 
Docp — 3w p a 2 r -4 / 3 ([Hansen et al.l 119751 IDaligaultl l2006f ) . where w 2 = ATrp(e/2m p ) 2 is the square of the ion plasma 

frequency, and a is the average ion spacing in the OCP. This yields Z?ocp = 0.03 cm 2 s _1 T^ 3 / p 31 ^ 54 , and thus nearly 
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the same result that diffusion is insignificant for t acc < 10 yr. 

REFERENCES 



Alastuey, A., & Jancovici, B. 1978, ApJ, 226, 1034 
Alcock, C, & Illarionov, A. 1980, ApJ, 235, 534 
Bahcall, J. N. 1964, ApJ, 139, 318 

Bildsten, L., Shen, K. J., Weinberg, N. N., & Nelemans, G. 2007, 
ApJ, 662, L95 

Bildsten, L., Townsley, D. M., Deloye, C. J., & Nelemans, G. 2006, 
ApJ, 640, 466 

Chapman, S., & Cowling, T. G. 1952, Mathematical Theory of 
Non-Uniform Gases (Cambridge; Cambridge University Press) 

Clayton, G. C, Geballe, T. R., Herwig, F., Fryer, C, & Asplund, 
M. 2007, ApJ, 662, 1220 

Daligault, J. 2006, PRL, 95, 065003 

Dearborn, D. S. P., Lattanzio, J. C, & Eggleton, P. P. 2006, ApJ, 
639, 405 

Eggleton, P. P. 1983, ApJ, 268, 368 

Fink, M., Hillcbrandt, W., & Ropke, F. K. 2007, A&A, 476, 1133 
Foley, R. J., et al. 2009, astro-ph/0902.2794 
Fujimoto, M. Y. 1982, ApJ, 257, 752 

Fujimoto, M. Y., Iben, I., Jr., & Hollowell, D. 1990, ApJ, 349, 580 
Fujimoto, M. Y., & Sugimoto, D. 1982, ApJ, 257, 291 
Gehrz, R. D., Truran, J. W., Williams, R. E., & Starrficld, S. 1998, 
PASP 110 3 

Glasner,'s. A., Livne, E., & Truran, J. W. 1997, ApJ, 475, 754 
Glasncr, S. A., Livne, E., & Truran, J. W. 2007, ApJ, 665, 1321 
Goriely, S., Jose, J., Hernanz, M., Rayet, M., & Arnould, M. 2002, 

A&A, 383, L27 
Gorres, J., et al. 1992, Nucl. Phys. A, 548, 414 
Graboskc, H. C, Dcwitt, H. E., Grossman, A. S., & Cooper, M. S. 

1973, ApJ, 181, 457 
Hansen, J. P., McDonald, I. R., & Pollock, E. L. 1975, Phys. Rev. 

A, 11, 1025 

Harris, M. J., Fowler, W. A., Caughlan, G. R., & Zimmerman, B. 

A. 1983, ARA&A, 21, 165 
Hashimoto, M., Nomoto, K., Arai, K., & Kaminisi, K. 1986, ApJ, 

307, 687. 
Iben, I., Jr. 1975, ApJ, 196, 525 

Iben, I., Jr., & Renzini, A. 1984, Phys. Rep., 105, 329 
Iben, I., Jr., & MacDonald, J. 1985, ApJ, 296, 540 
Iben, I., Jr., & Tutukov, A. V. 1989, ApJ, 342, 430 
Iben, I., Jr., & Tutukov, A. V. 1991, ApJ, 370, 615 
Iglesias, C. A., & Rogers, F. J. 1993, ApJ, 412, 752 
Iglesias, C. A., & Rogers, F. J. 1996, ApJ, 464, 943 
Iijima, T., & Nakanishi, H. 2008, A&A, 482, 865 
Itoh, N., Hayashi, H., Nishikawa, A., & Kohyama, Y. 1996, ApJS, 
102, 411 

Itoh, N., Totsuji, H., Ichimaru, S., & Dewitt, H. E. 1979, ApJ, 234, 
1079 

Karakas, A. I., Lee, H. Y., Lugaro, M., Gorres, J., & Wiescher, M. 

2008, ApJ, 676, 1254 
Kato, M., & Hachisu, I. 2003, ApJ, 598, L107 

Kato, M., Hachisu, I., Kiyota, S., & Saio, H. 2008, ApJ, 684, 1366 
Khokhlov, A. M. 1988, Ap&SS, 149, 91 
Khokhlov, A. M. 1989, MNRAS, 239, 785 
Kunz, R., et al. 2002, ApJ, 567, 643 

Lang, K. R. 1980, Astrophysical Formulae (Berlin : Springer) 
Limongi, M., & Tornambe, A. 1991, ApJ, 371, 317 



Livne, E. 1990, ApJ, 354, L53 
Livne, E. & Glasner, A. S. 1990, ApJ, 361, 244 
Livne, E. & Glasncr, A. S. 1991, ApJ, 370, 272 
Lodders, K. 2003, ApJ, 591, 1220 

Mocak, M., Miiller, E., Weiss, A., & Kifonidis, K. 2008, A&A, 490, 
265 

Mocak, M., Miiller, E., Weiss, A., & Kifonidis, K. 2008, astro- 

ph/0811.4083 
Naucnberg, M. 1972, ApJ, 175, 417 

Nelemans, G., Portcgics Zwart, S. F., Verbunt, F., & Yungelson, 

L. R. 2001, A&A, 368, 939 
Nomoto, K. 1982, ApJ, 253, 798 
Nomoto, K. 1982, ApJ, 257, 780 
Nomoto, K., & Sugimoto, D. 1977, PASJ, 29, 765 
Paquctte, C, Pelletier, C, Fontaine, G., & Michaud, G. 1986, 

ApJS, 61, 177 
Peterson, V. L., & Bahcall, J. N. 1963, ApJ, 138, 437 
Piersanti, L., Cassisi, S., Iben, I., Jr., & Tornambe, A. 1999, ApJ, 

521, L59 

Piersanti, L., Cassisi, S., & Tornambe, A. 2001, ApJ, 558, 916 

Piro, A. L., & Chang, P. 2008, ApJ, 678, 1158 

Prialnik, D., & Kovctz, A. 1984, ApJ, 281, 367 

Rogers, F. J., & Nayfonov, A. 2002, ApJ, 576, 1064 

Rosner, R., Alexakis, A., Young, Y.-N., Truran, J. W., & 

Hillebrandt, W. ApJ, 562, L177 
Saumon, D., Chabricr, G., & van Horn, H. M. 1995, ApJS, 99, 713 
Serenelli, A. M., & Fukugita, M. 2005, ApJ, 632, L33 
Serenelli, A., & Weiss, A. 2005, A&A, 442, 1041 

Shen, K. J., & Bildsten, L. 2007, ApJ, 660, 1444 

Shen, K. J., & Bildsten, L. 2009, ApJ, 692, 324 

Starrficld, S., Truran, J. W., Sparks, W. M., & Kutter, G. S. 1972, 

ApJ, 176, 169 
Sweigart, A. V., & Gross, P. G. 1978, ApJS, 36, 405 
Taam, R. E. 1980, ApJ, 237, 142 
Taam, R. E. 1980, ApJ, 242, 749 
Thomas, H.-C. 1967, Z. Astrophys., 67, 420 

Tilley, D. R., Weller, H. R., Cheves, C. M., & Chasteler, R. M. 

1995, Nucl. Phys. A, 595, 1 
Timmes, F. X. 1999, ApJS, 124, 241 
Timmcs, F. X., & Niemeyer, J. C. 2000, ApJ, 537, 993 
Timmes, F. X., & Swesty, F. D. 2000, ApJS, 126, 501 
Truran, J. W., & Livio, M. 1986, ApJ, 308, 721 
Tutukov, A., & Yungelson, L. 1996, MNRAS, 280, 1035 
van Teeseling, A., Reinsch, K., Hcssman, F. V., & Beuermann, K. 

1997, A&A, 323, L41 
Weinberg, N. N., Bildsten, L., & Schatz, H. 2006, ApJ, 639, 1018 
Woosley, S. E. 1986, in Nucleosynthesis and Chemical Evolution, 

ed. 6. Hauck, A. Maeder, & G. Magnet (Sauvcrny: Geneva 

Observatory) 

Woosley, S. E., Taam, R. E., & Weaver, T. A. 1986, ApJ, 301, 601 
Woosley, S. E., & Weaver, T. A. 1994, ApJ, 423, 371 
Yungelson, L. R. 2008, Astronomy Letters, 34, 620 
Yoon, S.-C, & Langcr, N. 2004, A&A, 419, 645 



